MAGNETIC PROPERTIES 
OF THE DOPED TWO-DIMENSIONAL ANTIFERROMAGNET 



A. SHERMAN 
Institute of Physics, Riia 142, 51014 Tartu, Estonia 

M. SCHREIBER 

Institut fur Physik, Technische Universitat, D-09107 Chemnitz 
and School of Engineering and Science, International University Bremen 
Campus Ring 1, D-28759 Bremen, Federal Republic of Germany 

The variety of the normal-state magnetic properties of cuprate high-T c super- 
conductors is interpreted based on the self-consistent solution of the self-energy 
equations for the two-dimensional t-J model. The observed variations of the spin 
correlation length with the hole concentration x, of the spin susceptibility with 
x and temperature T and the scaling of the static uniform susceptibility are well 
reproduced by the calculated results. The nonmonotonic temperature dependence 
of the Cu spin-lattice relaxation rate is connected with two competing tendencies 
in the low-frequency susceptibility: its temperature decrease due to the increasing 
spin gap and the growth of the susceptibility in this frequency region with the 
temperature broadening of the maximum in the susceptibility. 



Magnetic properties of cuprate perovskites have been extensively studied dur- 
ing the last years, both because of their unusual behavior and in the hope that 
they might provide insight into the physical origin of high-T c superconductivity.EjLl 
These properties provide a good test for models proposed for the description of 
cuprates and approximations used for their solution. One of the most popular 
models is the two-dimensional t-J model of Cu-0 planes. This model describes 
the low-energy dynamics of the realistic three-band Hubbard model.Q Recently we 
have suggested an approximate method for the solution of the t-J model, which has 
merits of retaining the rotation symmetry of spin componentsjn the paramagnetic 
state and of the absence of any predefined magnetic ordering.Q The obtained mag- 
netic properties of the model appear to be close to those observed in the cuprate 
perovskites which allows us to clarify mechanisms leading to the unusual behavior. 

The Hamiltonian of the two-dimensional t-J model reads 

H = ^ ] inmlij^miT + ~Z ^ ] ^nm ( s n s m s n s m ) — M ^ ' -^ni (1) 



n m 



where a na = |ner)(nO| is the hole annihilation operator, n and m label sites of 
the square lattice, a = ±1 is the spin projection, \na) and |n0) are site states 
corresponding to the absence and presence of a hole on the site. For nearest neighbor 
interactions t nm — t ^ a <5 n ,m+a and J nm = J J2 a ^n.m+a where t and J are the 
hopping and exchange constants and the four vectors a connect nearest neighbor 
sites. The spin-i operators can be written as s„ = \ cr|ner)(nc7| and s„ = 
|ner)(n, — oj, fi is the chemical potential and X n = |nO)(nO|. 

The method suggested in Ref. 7 is based on Mori's projection operator 
techniqueO which allows one to represent Green's functions in the form of con- 
tinued fractions and gives a way for calculating their elements. The residual term 
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of the fraction can be approximated by the decoupling which reduces this many- 
particle Green's function to a product of simpler functions. In this way we obtained 
the following self-energy equations for the hole G(k£) = — £0(i)({ak CT (t)> a k CT }) ana ^ 
spin D(kt) = — i#(i)([s|(i), sf_ k ]) Green's functions: 



D(kuj) = 



[4Ja(A + 1 + 7k )]- 1 n(kw) + 4JCi( 7k - 1) 



uj 2 - n(kcj) - 

G(kw) = 0[w — £k + /i — £(kw)]-\ 
where the self-energies II (kw) and S(kw) read 

ImII(kj) = — — (A + l + 7k )^( 7k -7 k+k 2 / <kj'[n F (u + Lj') -n F (co')} 

k' " / -°° 
Xi4(k + k^ W)A(kV), 

.(3) 



7 k-k' + 7k + sgn(w')( 7k -k' ~ 7k) \/ t — 

1 - 7k' 



x [iib(-w') + n F (w - u/)] A(k - k', a; - a/^kV). 

In the above formulas 7 k = \ X) a exp(ika), nj?(o;) = [exp(o;/T) + 1] — 1 , risiyj) = 
[exp(w/T) - l] -1 , A(koj) = -7r- 1 ImG(kw) and B(kw) = -7r -1 Im D(kw) are the 
hole and spin spectral functions, // = /i — 4.Fi0 -1 £ + 3G10" 1 J, = |(1 + a;) and 

^ = 16J 2 a|G 1 |(l- 7k )(A + l + 7k ), £k = (4# + 6G 1 ^ 1 i-3 J F 1 0- 1 J) 7k , (4) 

The hole concentration x, the spin Gi = (sj^sji) and hole i*i = (a|ai+ a ) 
nearest neighbor correlations are determined by the relations 

dujn F (uj)A(kuj), Gi = -^^7k y coth^— J ^(ku), 



Fi = ./V^^Tk / du;n F (uj)A(kcu). 

k 



(5) 



The parameter of vertex correction a is set equal to its value in. the undoped case, 
with introduced artificial broadening a = 1.802. As known,H in the considered 
two-dimensional system the long-range antiferromagnetic ordering is destroyed at 
any nonzero T and, as can be seen from the above formulas, at any nonzero x. 
The arising state with short-range order has to satisfy the constraint of zero site 
magnetization (sf ) = which can be written in the form 

V l - x ) = JjT, r^coth(^)s(M- (6) 

k J ° 

This equation fixes the parameter A describing the spin gap at (tt, n). 

The same derivation applied to the transversal spin Green's function ((sj^jsj^ 1 )) 
gives ((s k 1 js^ 1 )) = 2D(kuj) indicating that the approximation used retains properly 
the rotation symmetry of spin components in the paramagnetic state. 
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Figure 1. The dispersion of spin excitations. Vertical bars show decay widths |ImII(kaj)|/(2£J] t ). 
Points Y, M and S corresponds to k = (0, it), (it, it) and (tt/2, 7r/2), respectively. 



Figure 2. (a) Spin correlations vs. x for T = 0.02t. (b) Spin correlations along the diagonal [i.e., 
1 = (m, m)] of the crystal for x = 0.12. The respective temperatures are indicated near the curves. 



For low x and T the magnitude of the dispersion is small in comparison with 
t manifesting the band narrowing in the antiferromagnetic surrounding. 

Equations (||)-(||) form a closed set which can be solved by iteration. We carried 
out such calculations for a 20x20 lattice and the parameters t — 0.5 eV, J — 0.1 eV 
corresponding to cuprates. Small artificial broadenings — r\ and —2r\u>^ -q — 0.02i, 
were added to Im£(kc<j) and imll(kci;), respectively, to widen narrow lines and to 
stabilize the iteration procedure. 

The frequencies of spin excitations satisfy the equation 

lu 2 - Rcn(kw) - w£ = (7) 

[see Eq. (2)]. For low x and T their dispersion is close the dispersion of spin waves 
(see Fig. la.). The main difference is the spin gap at (tt, 7r) the magnitude of which 
grows with x and T (Fig. |l]b). In an infinite crystal this gap is directly connected 
with the spin correlation length £. Indeed, for large distances and low T we find 



5 i a o, 



= 7V- 1 ^e lkl |°°^coth(^)s(M oc (8) 



where Q = (jr, it) and £ = a/(2\/A) with the intersite distance a. For low x 
we found that A sa 0. 2x- and consequently £ w a/y/A. This relation has been 
experimentally observedE3 in La2_ 2; Sr 2 ;Cu04. 

As seen from Fig. [l]b, with growing x the spin excitation branch is destroyed in 
some region around the T point - Eq. (0) has no solution for real ui due to negative 
ReIT(kti;). For fixed T the size of this region grows with x. However, as follows from 
Fig. |l|b, with increasing T the branch is recovered in this region. This result is 
consequence of the temperature broadening of the spin-polaron peak in A(kcu) 
Large negative values of ReIT(ka;) which are the reason for the lack of a solution of 
Eq. (^) are connected with this peak. With its broadening |ReII| becomes smaller. 

The spin correlations C m i = (sfsg), 1 = (n,m), Eq. (||), are shown in Fig. ||. 
Distant correlations decrease rapidly with x. For large enough x and T the corre- 
lations decay exponentially with distance in the considered finite lattice. 

The magnetic susceptibility is connected with the spin Green's function (^|) by 
the relation x z (kw) = — 4/i^D(kaj), where \jlb is the Bohr magneton. Experiments 
on inelastic neutron scattering give information on the dynamic spin susceptibility 
which can be directly compared with the calculated results. Such comparison is 
carried out in Fig. [| YBa2Cu30y_j, is a bilayer crystal and the symmetry allows 
one to divide the susceptibility into odd and even parts. For the antiferromagnetic 
intrabilayer coupling the odd part of the susceptibility can be compared with the 
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Figure 3. The imaginary part of the spin susceptibility for k = (it, it). Curves demonstrate results 
of our calculations for T = 0.024, and x = 0.043, 0.08 and 0.12 (from top to lptrttom). Squares 
show the imaginary part of the odd spin susceptibility for k = (it, it) measureocl in normal state 
YBa 2 Cu307_ H at T = 100 K for y = 0.5, 0.17 and 0.03 (from top to bottom). 



Figure 4. The uniform static spin susceptibility vs. T (a) and x (b). 



calculated results. The oxygen deficiencies y = 0.5, 0.17 and 0.03 in the experi- 
mental data-in Fig. || correspond to the hole concentrations x = 0.07, 0.13 and 0.17 
respectively.til We scale the experimental Im^ so that their maximal values coin- 
cide approximately with the calculated maxima. In absolute units the calculated 
maximal values are approximately 2-2.5 times larger than the experimental ones. 
As seen from Fig. ^ the calculated data reproduce correctly the frequency depen- 
dence of the susceptibility, the values of the frequency for which Imx(7T, 7r) reaches 
maximum and their evolution with doping. We notice also that the calculated 
temperature variation of the susceptibility is in good agreement with experiment. 

The calculated dependencies of the uniform static spin susceptibility Xo — 
x(k — > 0, a; = 0) on T and x are shown in Fig. [|. In this figure xo/Vn-|hcs 
in the range 2-2.6 eV™ 1 which is close to the values 1.9-2.6 eV _1 obtaincda for 
YBa2Cu307_ a . The temperature dependence of xo has a maximum and its tem- 
per ature T m grows with decreasing x. Analogous behavior of xo was observedBB in. 
cuprates for large enough x. In Fig. f|a T m w 600 K which is close to T m observed^ 
in La2-xSr x Cu04 for comparable x. As known, in the undoped antiferromagnet 
T m w J. On the high-temperature side from the maximum Xo(T) tends to the 
Curie- Weiss dependence 1/T. As seen from Fig. |]a, the two curves for the different 
x are very close in shape and can be superposed by scaling to the same maxi- 
mal value of x& an d T m . Analogous scaling was observedEJ in La 2 _ 2; Sr :I .Cu04. As 
in experiment^ the dependence Xo( x ) m Fig- @b has a maximum. However, the 
experimental maximum is achieved at x ~ 0.25. 

It is of interest to clarify the reason for the nonmonotonic behavior of the de- 
pendence xo(T). From Eqs. (||) and (||) it is clear that x(k,w = 0) is continuous 
for small k. Therefore let us consider a small but finite value of k for which 
X(k,0) cx J™ du/BQuSyw'. The maximum in B{ku') is shifted to lower fre- 
quencies and its intensity decreases with increasing temperature. However this 
dependence is superimposed with the function 1/lo' in the integral. Nonmonotonic 
behavior of Xo(T) is the consequence of the existence of some optimal location of 
the maximum in B(k.Lu') when the area under the curve B(ku)')/u>' is maximal. 
Johnston's scalingp means that the temperature softening of the long-wavelength 
maximum occurs similarly for different x and that holes and temperature fluctua- 
tions act in a similar manner as sources of the softening. 

The spin-lattice relaxation and spin-echo decay rates were calculated with the 
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Figure 5. The temperature dependencies of the spin-lattice relaxation and spin-echo decay rates. 
Open circles with right axes represent experimental results, filled circles ™jth left axes are our 
calculations^ (a,c,d) calculations for H||c and x = 0.12, measurements!] in YBa2Cu30n.63 
(x « 0.1).Lj (b) calculations for nonoriented configuration with x = 0.043, measurementsn in 
Lai.9 6 Sr .04CuO4. 



use of the equation 
1 



a T lfs T 2n 2 B N 
1 0.69 



63T| G 128/4 I N k 



W , w -0, 

(9) 

J= £ 63 F e 2 (k) [Re X (kO)] 2 - 1^ 63 ^(k)Re X (kO) 



AT 

k 



where the form factors "^(k) and hypcrfine coupling constants were taken from 
Ref. 4. The superscript a = 63 or 17 indicates that the respective quantity belongs 
to Cu and O, respectively. The subscripts f3 = || or _L refer to the direction of the 
applied static magnetic field H with respect to the axis c perpendicular to the Cu-0 
plane. The form factor 63 F e is the filter for the Cu spin-echo decay time 63 T2G- Our 
calculated results and the respective experimental data are given in Fig. g. The 
calculations reproduce satisfactorily main peculiarities of the temperature depen- 
dence of the spin-lattice and spin-echo decay rates. The growth of (^TiT) -1 with 
decreasing x is connected with the increase of the spectral intensity of spin exci- 
tations near (tt,tt). For the same x ( 63 TiT) _1 is one to two orders of magnitude 
larger than ( 17 TiT) _1 . This is a consequence of the spin susceptibility which is 
strongly peaked near (tt, tt) and the form factors which test different k regions.Q 

In Fig. H the calculated values of the spin-lattice relaxation rates are nearly an 
order of magnitude smaller than the experimental ones. This difference is connected 
with the approximation made in the calculation of D(ku>) the continued fraction of 
which was cut short in the second link. This leads to an additional power of ui in 
the denominator of B(kuj). 

It is of interest to clarify the reason for the nonmonotonic temperature de- 
pendence of ( 63 TiT) _1 for moderate x (see Fig. ||a). The main contribution to 
( 63 TiT) -1 is given by the vicinity of the (tt, tt) point. For low oj in this region 
Imx(kw) first grows then decreases with increasing T. This is connected with two 
competing tendencies in the spin spectral function: i) the spin gap A grows with 
T which leads to a decrease of B in the mentioned regions of k and u> and ii) the 
temperature broadening of the maximum in B which leads to its growth in these 
regions. For moderate x and low T practically A does not depend on T, the second 
tendency dominates and ( 63 TiT)~ 1 grows with T . For higher T the gap grows more 
rapidly with T leading to the predominance of the first tendency and ( 63 TiT) _1 
decreases with T, For small x temperature fluctuations dominate, A grows notice- 
ably even for small T and ( 63 TiT) -1 decreases monotonously with T (see Fig. |^b). 
Due to the form factor 17 F||(k), the central part of the Brillouin zone makes the 
main contribution to ( 17 TiT) _1 . For moderate T and x the existence of the spin 
gap affects B(ku) weakly here. For small li the temperature dependence of B(kuj) 
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is determined by broadening. As a consequence, ( 17 TiT) _1 grows monotonously 
with T (see Fig. ||c). The low-temperature decreases of ( 63 TiT) _1 and x° were 
considered as manifestations of the spin gap. As follows from the above discussion, 
these peculiarities are not directly connected with the gap. 

This work was partially supported by the ESF grant No. 4022 and by DFG 
(SFB 393). 
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